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ABSTRACT 

We present a theoretical delay time distribution (DTD) of Type la supernovae on the basis of our 
new evolutionary models of single degenerate (SD) progenitor systems. Our model DTD has almost 
a featureless power law shape (oc with n « 1) for the delay time from t ^ 0.1 to 10 Gyr. This 
is in good agreement with the recent direct measurement of DTD. The observed featureless property 
of the DTD has been suggested to be favorable for the double degenerate (DD) scenario but not for 
the SD scenario. If the mass range of the companion star to the white dwarf (WD) were too narrow 
in the SD model, its DTD would be too limited around the companion's main-sequence lifetime to 
be consistent with the observed DTD. However, this is not the case in our SD model that consists of 
the two channels of WD -I- RG (red giant) and WD + MS (main-sequence star). In these channels, 
the companion stars have a mass range of ^ 0.9 - 3 Af© (WD+RG) and ^ 2 - 6 Af© (WD-I-MS). 
The combined mass range is wide enough to yield the featureless DTD. We emphasize that the SD 
scenario should include two important processes: the optically thick winds from the mass-accreting 
WD and the mass-stripping from the companion star by the WD wind. 

Subject headings: binaries: close — galaxies: evolution — stars: winds, outflows — supernovae: 
general 



1. INTRODUCTION 

Type la supernovae (SNe la) play the important roles 
in astrophysics as a standard candle to measure cos- 
mological distances as well as the production site of 
a large part of iron group elements. However, the 
nature o f SN la progenitors has not been clarified 
yet ( e .g..lHillebrandt fc Niemever"! 120001: iNomoto et al. I 
119971 l2000HLivioll200Gf ). It has been commonly agreed 
that the exploding star is a carbon-oxygen white dwarf 
(C-l-0 WD) and the observed features of SNe la are bet- 
ter explained by the Chandrasekhar mass model than 
the sub- Chandrasekhar mass model. However, there has 
been no clear observational indication as to how the WD 
mass gets close enough to the Chandrasekhar mass for 
carbon ignition; i.e., whether the WD accretes H/He- 
rich matter from its binary companion [single degenerate 
(SD) scenario] or two C-l-0 WDs merge [double degen- 
erate (DD) scenario]. 

It has been suggested that SNe la have a wide range 
of delay time from t < 0.1 Gyr to t > 10 Gyr (e.g., 
iMannucci et al. I I2006D . Here the delay time, t, is de- 
fined as the age at the expl osion of the SN la proge nitor 
from its birth. According to IMannucci et alH (|2006f ). the 
present observational data of SNe la are best matched by 



a bimodal population of the progenitors, in which about 
50 percent of SNe la explode soon after their stellar birth 
at the delay time of t ~ 0.1 Gyr, while the remaining 50 
percent have a much wider distribution of the delay time 
around i ~ 3 Gyr. If the delay time distribution (DTD) 
of SNe la is observationally obtained, we are able to pre- 
clude some models that are inconsistent with the DTD. 

Recently, the direct measureme nt of the DTD has been 
reported bv iTotani et al. I (|2008[ ). Their DTD shows a 
featureless power law (oc n 1) between t ~ 0.1 
and ~ 10 Gyr. On the basis of their results, they 
argued that the DTD strongly supports the DD sce- 
nario of SN la progenitors mainly because the feature- 
less power law distribution is in good agreement with 
the prediction of the DD scenario. They also concluded 
that the SD scenario is not well supported mainly be- 
cause some "detailed" binary population synthesis codes 
predict prominent peaks in the DTD at characteristi c 
time scales (e.g.. iBelczvnski et al . 1 120 051: iMeng et al" 
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20081 iRuiz-Lapuente fc Canal 1119 98: Yu ngelson fc Livio 
2000^. although some "simple" SD models with sim- 
plified treatments of binary evolution ha ve broad DTD 
shapes similar to the observed DTD fe.g..[Greggio ]|2005l : 
iMatteucci et al. ll2006HKobavashi fc Nomoto \\20m . 

In the SD scenario, the DTD is closely related to the 
main-sequence (MS) lifetime of the companion star and 
thus the initial mass of the companion (secondary), M2fl- 
This is because the mass transfer from the companion to 
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the WD starts when the companion evolves off the MS 
and expands to fill its Roche lobe. The mass-accreting 
WD quickly grow s to the Chandrasekha r mass limit in 
~ 10^ yr (see e.g.. iHachisu et al. Ill999allbl . l2008f) . which 
is much shorter than the age of the binary system t > 
10* yr. The featureless power law (oc t~") of the DTD 
between t ^ 0.1 and ^ 10 Gyr requires a very broad 
mass distribution of the companion from M2 ^ 0.9 to 
~ 6 Mq. 

In the "detailed" SD models, they follow each binary 
evolution including many binary evolutionary processes 
and, as a result, the mass of the companion that can 
produce an SN la is constr ained to a certai n range, e.g., 
M2,o « 2 - 3.5 Mq in Me ng et al. I (|2008D . Therefore, 
the delay time is also constrained to a narrow range of 
t ~ 0.3 — 1.2 Gyr, which is inconsistent with the observed 
DTD. 

In some "simple" SD models, a mass range of the com- 
panion is assumed a priori to be M2,o = 0.8 — 8 AIq 
simply from the condition of M2,o < -^1,0 without tak- 
ing into account the constraints on the mass accretion 
rate onto the WD and thus on the separat ion between 
the WD and the companion (jGreggio 1 120051 ). 

In the present paper, we show that the SD model 
with taking into account the "detailed" binary evo- 
lution is viable against the observed DTD. Actually, 
the required broad distribution of the companion mass, 
M2,o ^ 0.9—6 Mq has already been predi cted by a recent 
"deta iled" progenitor model of SNe la (jHachisu et al. I 
I2OO8I ). However, Hachisu et al. did not present any 
DTDs mainly because at that time there were no ob- 
servational data to compare with the theoretical results. 
Here we pre sent DTDs on the bas is of the new "detailed" 
SD models (|Hachisu et al. I2008D in ^^2 and 3 and com- 
pare with the observation ( Totani et al. |[2008[ ) in §4. 
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2. MASS-STRIPPING EFFECT AND BINARY EVOLUTION 

We start the binary evolution from the zero-age MS. 
Unless the initial separation of the binary components 
is too close, the more massive (primary with the mass 
of Afi^o) component evolves to a red giant star (with a 
helium core) or an AGB star (with a C-l-0 core) and fills 
its Roche lobe. Subsequent mass transfer from the pri- 
mary to the secondary is rapid enough to form a common 
envelope. The binary separation shrinks greatly owing 
to the mass and angular momentum losses from the bi- 
nary system during the first common envelope evolution. 
The hydrogen-rich envelope of the primary component is 
stripped away and the primary becomes a helium star or 
a C-l-0 WD. The helium star further evolves to a C-l-0 
WD after a large part of helium is exhausted by core- 
helium-burning. Thus we have a binary pair of the C-l-0 
WD and the secondary star that is an MS star; the mass 
of the secondary star, M2, is still close to ^^2,0, because 
the accreted mass during the common envelope phase is 
negligibly small. 

After the secondary evolves to fill its Roche lobe, the 
WD accretes mass from the secondary and grows to 
the critical mass (Mja — 1.38 Mq) and explodes as 
an SN la if the initial binary orbital period (Pq) and 
the initial mass of the secondary (M2,o) are inside the 
regions (labeled "initial") shown in Figure [TJ There 
are two separate regions; one is for binaries consist- 
ing of a white dwarf and a main-sequence star (WD 
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Fig. 1. — The regions that produce SNe la are plotted in the 
log P — M2 (orbital period - secondary mass) plane for the white 
dwarf and main-sequence star (WD + MS) system (left) and the 
white dwarf and red giant (WD + RG) system (right). Here we 
assume the metallicity of Z = 0.02 and the initial white dwarf 
mass of A'/vvD = t-0 The initial system inside the re- 

gion encircled by a thin solid line (labeled "initial") increases its 
WD mass up to the critical mass (Mia = 1.38Mq) for the SN 
la explosion, the regions of which are encircled by a thick solid 
line (labeled "final"). Currently known positions of the recur- 
rent novae and supersoft X-ray sources are in dicated by a star 
mark (*) for U Sco (e.g.. ^ Hachisu et al. 11200031 ). a triangle for T 
CrB (e.g. [Be lczynski & Mikolaicwska 1998'), a square for V Sge 
JHachisu fcKato 2003 c). but by arrows for the other three recur- 
rent novae, V394 CrA, CI Aql, and RS Oph, because the mass of 
the companion is not yet available explicitly. Two subclasses of the 
recurrent novae, the U Sco type and the RS Oph type, correspond 
to the WD + MS channel and the WD + RG channel of SNe la, 
respectively. 

-I- MS) and the other is binaries consisting of a white 
dwarf and a red giant (WD -I- RG). Here the metallic- 
ity and the initial white dwarf mass are assumed to be 
Z = 0.02 and AfwD,o = 1-0 Mq. Note that the WD with 
MwBfi — 1-0 Mq for Z = 0.02 forms from the primary 
star of Mi,o 7 Mq (jUmeda et al. Ill999f ). 

Here we emphasize the effects of two important pro- 
cesses in the binary evolutions. The first one is the 
accretion wind evolution. Mass-accreting WDs blow 
strong winds when the mass transfer rate onto the WD 
exceeds the critical r ate of Mcr 1 x IO^^Mq yr"-'^ 
(jHachisu et al. I [1996? ). The angular momentum taken 
away by this fast wind is much smaller than the orbital 
angular momentum, so that the separation of the binary 
does not shrink (jHachisu et al. |[T999a ) . If we ignore this 
wind effect, almost all the binaries would undergo a com- 
mon envelope phase and merge (the WD -I- MS systems) 
or form a double degenerate system (the WD -I- RG sys- 
tems) due to a large amount of angular momentum loss. 
None of them becomes an S N la through the SD channel 
(e.g.. IHachisu et al. |[r999a[ ). 

The second one is the mass-stripping from the sec- 
ondary surface by the WD winds. This attenuates the 
mass transfer rate from the secondary to the WD, so that 
the binary can avoid the formation of a common enve- 
lope even for a rather massive secondary (jHachisu et al. I 
Il999al I2OO8D . In our results in Figure [U this mass- 
stripping effect is critically important. 

Our results showed that the "initial" region of WD 
-I- MS systems extends up to such a massive (M2,o ~ 
5 — 6 Mq) secondary, which consists of a very young 



Delay Time Distribution 



3 



>- -1 



' ' (o) 'f(q)=1 



WD+MS systen 




WD+R 



-1 -0.5 0.5 

log (delay time) (Gyr) 



» -1 



(b) f(q) = 2/(1+q)2 



WD+MS system/ 




-0.5 0.5 1 

og (delay time) (Gyr) 



(a) f(q) = 24qV(1+q)' 




WD+MS system 



WD+R 



sysl em 



-0.5 







0.5 



log (delay time) (Gyr) 



(d) a=2.35, f(q) = 1 




WD+MS systen 



WD+R 



-0.5 







0.5 



og (delay time) (Gyr) 



Fig. 2. — The delay time distributions (DTDs) for our single degenerate (SD) models of SNe la. The ordinate is the DTD in units of per 
century and per 10^^ Lk.P)- Open circles with an open box are the observational DTD taken from lTotani et al. I fOOOB ) (for < 10 Gyr) and 
IMannucci et al.~l l |2005f ) (at 11 Gyr). Each open box indicates a one-sigma error of each measurement. Throughout all figures from (a) to 
(d), thick solid lines indicate the total DTD coming both from the WD + MS and WD + RG systems. Each contribution from each channel 
is separately shown by thin solid lines. The data of our DTD only for (a) are tabulated in Table [T] The binary mass ratio distribution and 
the initial mass function power are assumed to be (a) f{q) = 1 and a = 2.5, (b) /(g) = 2/(1 + g)^ and a = 2.5, where a factor of 2 is a 
normalization factor for Jq f{q)dq = 1, and fwD+MS = 0.0038 yr"^ and u^-d+rg = 0.0038 yr~^ in our Galaxy, (c) f{q) = 245^/(1 + g)* 
HGregeio &: Renzini 11198^) and a = 2.5, where a factor of 24 is a normalization factor for Jq f(q)dq = 1, and !AvD+MS = 0.0034 yr~^ and 
^'WD+RG = 0.0025 yr~^ in our Galaxy, and (d) f{q) = 1 and a = 2.35. Dotted line: same as that in Fig. [2^, dashed line: same as that in 
Fig. [2}3, dash-dotted line: same as that in Fig. [5];. 



population of SNe la with such a short delay time as 
t<0.1 Gyr. On the other hand, the WD + RG systems 
with a less massive RG (M2,o ~ 0.9 — 1.0 Mq) consist of 
a very old population of SNe la of < > 10 Gyr. 

3. DELAY TIME DISTRIBUTION 
The birth rate of SNe la in our Galaxy is expressed as 
dMi,o 



^ = 0.2 



^2.5 



f{q) dq dlog a yr" 



(1) 



ID {Mi,ol 

along equation (1) of llben fc Tutukov I (|1984[ ). where 
q = M2,o/Mifi, a, and Mi^ are the mass ratio, the 
separation, and the primary mass in solar mass units, 
respectively, at the birth of a binary and D is the SN 
la region in the (Mi_0: 9, a)-space. Dividing the three- 
dimensional space of (Mi_o, q, a) into grids of {Mi, qj, ak), 
we follow the binary evolution starting from the initial 
state of {Mi,qj,ak) and obtain v for the binaries which 
finally explode as SNe la. At the same time, we record 
the delay time of t for each {Mi, qj, ak)- In our calcula- 
tion, the SN la region of D consists of these SN la grid 



points. This procedure can easily be done when we know 
the "initial" SN la regions for both the WD + MS and 
WD -I- RG systems as shown in Figure [T] 

In our prev i ous works as seen in Figure 4 of 
iHachisu et aTl (I2008D for t he WD + MS system and 
Figure 12 of iHachisu et all (fl999al) for the WD + RG 
system, we have obtained such SN la regions for differ- 
ent initial WD masses, Mwd,o = 0.7, 0.8, 0.9, 1.0, and 
1.1 Mq, which formed from the primary star of masses 
Mi^o ^ 4, 5, 6, 7, and 8 Mq at t he birth, respectively, 
for'z = 0.02 (jUmeda et al. II1999D . Using these resuhs, 
we have estimated the SN la birth rate in our Galaxy 
as i^wD+MS = 0. 0035 yr~^ and y wu+HG = 0.0032 yr~^, 
respectively see IHachisu et aT~i ri999a.b. 200 8), for more 
detail). Here we assume f{q) — 1 along I ben fc Tutukov I 
(|l984f ). Therefore, the relative ratio of the SNe from the 
WD -I- MS progenitors to those from the WD + RG pro- 
genitors is tms/rg = l-l- 

Now we estimate the DTD of SNe la forming from 
the WD + MS and WD + RG systems, by integrating 



Hachisu et al. 



4 

only the initial sets of (M^, qj, a^) having the delay time 
between t — and t + At, as 

(2) 

where we adopt f{q) = 1 and 10 bins of the delay time 
at t = 0.05, 0.1, 0.2, 0.4, 0.8, 1.6, 3.2, 6.4, 12.8, and 25.6 
Gyr. The resultant DTD as summarized in Table [U is 
not normalized, thus being scale-free. We normalize our 
DTD by fitting the value to the observation at 11 Gyr 
as shown in Figure [2^. Our DTD shows a featureless 
power law (c>c t~", n k, 1) from 0.1 to 12 Gyr, which is 
in good agreement with Totani et al.'s (2 08) D TD and 
the data at 11 Gyr bv iMannucci et al. I (|2005() . Here, 
we assume = 70 km s~^Mpc~^ for the Hubble con- 
stant. Thus the DTD on the basis of our "detailed" SD 
model of SN la progenitors is consistent with the obser- 
vation. It should be noticed that we plotted the number 
ratio of SNe la aga inst the delay time in Figure 12 of 
iHachisu et al. I (|2008i) . which clearly shows a bimodality 
of SN la progenitors, but the number ratio itself is not 
equal to the DTD of SNe la because the DTD is time- 
derivative of the number ratio. 

Here, for instructive purposes, we derive an approx- 
imate power law of our DTD. The main-sequence life- 
time of the secondary, i2, can be estimated as t2 oc 
M2,o/L2,o fx {M2,o)^-"\ where the mass-luminosity rela- 
tion at the zero-age main-sequence is approximately writ- 
ten as 2.2,0 oc (Af2,o)" with m = 3.5 for M2,o = 3-7 Mq. 
The appropriate range of initial separation (A log a) is 
roughly proportional to the range of orbital period, i.e., 
A log a « (2/3) A log P, so the area of AM2,o A log a is 
calculated approximately from the SN la region in the 
logP — M2 plane for a given Mi.o like in Figure [T] This 
area becomes narrower as tli6 initicil WD mass, -/ V^wp p, 
decreases like Figure 4 of IHachisu et al. I (|2008[ ). The 
initial WD mass is closely related to the initial primary 
mass, Mi,o, so that we numerically obtain the approxi- 
mate relation of 




dM2,o dlog a oc (Ml,o)'■^ (3) 



for the fixed Mi^ between 6 < Mi^ < 9. We also nu- 
merically obtain a similar relation of 




(Mi,o)-^-'dMi,o dlog a cx (A/2,o)^'-°, (4) 



for the fixed A'hfi between 3 < M2,o < 6. Then the DTD 
can be approximated as 

2AtJt_^t M2,Q 2AtJ,_^, t2 

for At ^ ^2- Here, we assume f{q) — 1 and Af2,o ^ 
(h)^^^^'^- Strictly speaking, the power of M2,o in equa- 
tion ^ is somewhere between —1.0 and —0.5. If we 
adopt the power of —0.5, the final power of the delay 
time, t, in equation ^ changes from —1.0 to —1.2 

For the WD -|- RG system, we use m — 5 for the mass- 
luminosity relation of the zero-age main-sequence stars 



with M2 = 0.7-2 Mq and Afa.o oc (12)^'^/^. Applying 
the area of the SN la regions th a.t have been already ca l- 
culated as shown in Figure 12 of lHachisu et al. I ()1999aD . 
we numerically obtain the same approximate relation as 
given by equation Q for 7 < Afi,o < 9 and equation ([4|) 
for 0.9 < Af2,o < 2. Then we have the same power law 
index as in equation ([5]), i.e., DTD(t) oc for the WD 
-I- RG channel. 

In both the WD + MS and WD + RG channels, the 
DTD has the power law index (oc t~") close to n = 1 
regardless of the mass-lifetime dependence of M2,o oc 

(^2)"^/^-^ for the WD + MS or M2,o oc (^2)"^^'' for the 
WD -I- RG system. The important "details" to realize 
such DTDs is how the SN la region shrinks or expands 
as the initial secondary mass, Af2,o, decreases. Here, the 
approximate relation given by equation ([4]) holds for both 
the WD + MS and WD + RG systems, which leads to 
the logarithmic form of the DTD (t) oc J d\ogM2fi/ At oc 
J d\ogt2/ At oc t~^ regardless of the mass-lifetime de- 
pendence, as seen in equation ([5|). 

4. CONCLUSIONS AND DISCUSSION 

iTotani et afl (|2008D argued that the "detailed" SD 
models on the basis of the population synthesis should 
have prominent peaks in the DTD at characteristic time 
scales of the secondary mass, thus being inconsistent with 
the observation. As already shown in Figure[2^, however, 
our DTD on the basis of the "detailed" binary evolution 
models has a featureless power law, being in good agree- 
ment with the observation. This is because the mass 
of the secondary star of the SN la system ranges from 
M2.0 ~ 0.9 to 6 Mq (Fig. [J) due to the effects of the WD 
winds and the mass stripping. In our model, moreover, 
the number ratio of SNe la between the WD -I- MS com- 
ponent and the WD -I- RG component is Tms/rg = 1.1- 
Such almost equal contributions of the two components 
help to yield a featureless power law as discussed below. 

As for the met allicity effect (e.g., 
iKobavashi &: Nomoto I l2008f l. we assume that the 
metallicity had already increased to Z = 0.02 (or more) 
at the birth of prog enitor stars mainly because the 
galaxies studied by , Totani et al. I ([2008) consist of old 
galaxies, the metallicity of which had already increased 
to Z = 0.02 or more. Therefore, Totani et al.'s (2008) 
data should not show any metallicity effect even if the 
metallicity effect really exists. 

In order to see the effect of different initial distributions 
of the binary mass r atio q, we calculate the DTDs for two 
cases of f(q) (e.g.. iGreggio &: Renzini I Il983l : iGreggio I 
2005; Kobavashi fc Nomoto II2008D . Figureia^ shows that 
the DTD for /(g) = 2/(1 -I- g)^ is in good agreement with 
the observation. The other DTD for f{q) = 24g2/(i + g)4 
shown in Figure [2j; is marginally consistent with the ob- 
servation. 

Such a weak dependence stems from the fact that the 
ratio of the two SN la components, ^ms/rGi depends 
slightly on f{q). The WD -I- MS component has a 
relatively short delay time, and tends to have a large 
M2,o- In contrast, the WD -|- RG component has a 
long delay time, having a small M2,o- In Figure [U 
for example, Mi.o ^ 7 Mq for Mwd.o ^ 1 Mq, while 
M2,o ^ 0.9 - 3 Mq for the WD -I- RG system, i.e., 
g- 0.13-0.4. The g distribution of /(g) = 24gV(l + g)4 
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TABLE 1 
Delay time distribution of SNe Ia'' 



delay time 
(Gyr) 



WD + MS'' WD + RG" 



total 



in our Galaxy 

has a peak at g = 1 and takes smaller values at smaller q. 
As a result, tms/rg = 1-4 for this f{q) = 2^"^ / {1 + q)^ is 
slightly larger than rMs/RG = 1-0 for f{q) = 2/(1 + g)^. 
Thus the DTDs of our "detailed" SD models is not so 
sensitive to the mass ratio distribution f{q), suggesting 
that almost a featureless power law shape (oc t~'^, n « 1) 
from i ~ 0.1 to 10 Gyr is common among our "detailed" 
SD models as long as the mass range of the secondary is 
M2,o ~ 0.9 - 6 Mq. 

We also calculate a DTD with a different power law 
index of IMF, i.e., a — 2.35 instead of a = 2.5 in equa- 
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tions (III) and As shown in Figure the results 
are hardly affected by changing a of the IMF as long as 
a = 2.35 ±0.15. 
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